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^ ■ Abstract 



^ ■ I discuss a formahsm for computing quantum scattering amphtudes using 

. a semiclassical expansion of a functional integral representation for the S- 

matrix. The classical background for the expansion is determined by solving 
\Q • the equations of motion subject to nontrivial boundary conditions deter- 

O . mined by the initial and final quantum states. The formalism is designed 

to accomodate intrinsically nonperturbative processes such as baryon num- 
ber violation, quantum tunneling and multiparticle scattering. It appears 
to yield a controlled small coupling expansion even at asymptotically high 
energies where instanton methods fail. 
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: INTRODUCTION 

X 

d • The behavior of quantum scattering amplitudes at large center of mass energies, 

E E^, = m/ g"^, is an open question currently under investigation. Here, m is a mass 
scale characterizing the particle excitations of the theory, and g is a. small coupling con- 
stant controlling the nonlinear interactions. It has been speculated that multiparticle 
scattering from initial states with only a small number of particles could be unsup- 
pressed at these energies. A second related issue is whether anomalous processes, such 
as baryon number violation in electroweak theory where E^ c:^ O {Mw/aw) ~ 10 TeV , 
are unsuppressed at energies of the same order |l[]. 

Standard perturbation theory is clearly inadequate to address either of these issues. 
Recent investigations have focused on expansions about nonperturbative solutions of 
classical field equations, such as instantons |]1], |^. However, none of the previously 
proposed solutions provides a stationary point for a controlled semiclassical expansion 
at the energies of interest, ~ i?* [Q. In general, the semiclassical expansion about 
the proposed instanton-like saddlepoints is controlled only at low energies and provides 
little or no information at E' ~ i?*. 

In this talk, I will discuss a method for finding classical solutions which will prove 
suitable saddlepoints of scattering amplitudes even at ~ E,,. The equations and 
boundary conditions which must be satisfied are dramatically simplified compared to 



*Talk presented at Sintra '94 workshop on Electroweak Physics and the Early Universe. 



a previous approach to the same problem by Mattis, McLerran and Yaffe [Q. In 
particular, the resulting equations and boundary conditions for the classical saddlepoint 
are no longer of the integro-differential type. Rather, the equations of motion are 
source-free. 



and in certain cases the boundary conditions are of the initial value type, which can 
be integrated forward in time using standard numerical algorithms. This is also a 
departure from the approach of Khlebnikov, Rubakov and Tinyakov |]^, ^ , who derive 
similar source- free equations, but for an initial coherent state, rather than a wave-packet 
type state which is relevant for high energy collisions. 

In light of this new formalism, I will discuss the implications of some recent nu- 
merical simulations in classical field theory for high energy two-particle scattering. In 
particular, in 4-dimensional pure SU(2) Yang-Mills theory. Gong, Matinyan, Miiller and 
Trayanov |^ have found an instability of high frequency standing waves of pure glue 
to decay to low frequency modes. When interpreted within our formalism, this result 
leads to the speculation that multiparticle scattering in SU{2) gauge theories may be 
unsuppressed (at least exponentially) at high energies. A generalization of our method 
to complex trajectories which allows their application to essentially any wave packet 
- coherent state scattering amplitude. I finally discuss how our formulation relates 
to some previous work by other authors on high energy multiparticle and anomalous 
processes in electroweak theory p, ^. 

Before launching into the formalism, let me make a few remarks on the method 
and the results we can expect. We are interested in quantum scattering amplitudes 
(or S-matrix elements) between arbitrary initial and final states \i) and |/). In most 
of what follows in the talk I will be specifically interested in states that are relevant to 
accelerator experiments (i.e. wave packet states), but here I will be completely general. 

Now, imagine thinking about the scattering process in terms of paths inside the 
functional integral. The quantum fields do whatever they want, weighted by the appro- 
priate action exp[iS], but also by the overlap of the asymptotic part of the path with 
the initial and final states. Therefore we expect that the amplitude must be expressible 
in the form: 

~ j d<pa<Pf D<P^,[<p{Ti)]^f[(t){Tf)]e''^^\ (1) 

where I have explicitly written the measure of path integration to include fluctuations 
of the fields in the asymptotic past and future. The wave-functionals '^ij measure the 
overlap of the initial and final states with position (or field operator) eigenstates </> at 
asymptotic times. Below I will give an actual derivation of (|l]) along with the explicit 
form of the wavefunctionals in the case that the initial and final states are either wave- 
packets or coherent states. (Recall that a plane wave is a special type of wave packet 
that is narrowly peaked in momentum space, so we can work with plane waves as well 
if we wish.) 

Now, the usual way to compute (|I]) at weak coupling is by expanding the functional 
integral about either the vacuum - which generates standard perturbation theory - or 
about a more exotic classical path such as the instanton, but which also happens to 
obey vacuum boundary conditions. (In the case of the instanton ([^) is first continued to 
Euclidean space, where finite action solutions must satisfy vacuum boundary conditions. 
The results of expanding the path integral are then analytically continued back term 
by term. More on this later.) Clearly, for some initial and final states it is not a big 
loss to evaluate the overlaps ^E'j,/[0] on the asymptotic vacuum. We know that in most 
instances standard perturbation theory works quite well. We also know that instanton 
methods give a controlled way of computing tunneling phenomena - at least at low 
energies. 



However, to get a better approximation to the above path integral, one could try 
to extremize the entire integrand, which will yield some nontrivial boundary conditions 
on the classical solutions, as determined by the specific form of the wavefunctionals, or 
in other words the initial and final states of interest. 

FORMALISM 



As I mentioned above, we wish to calculate the S-matrix element between an 
initial two particle state and an arbitrary final state tO. The result will be a boundary 
value problem which determines the classical field which is the stationary point of the 
S-matrix element. 

To this end, we express the kernel of the S-matrix in a basis of coherent states, first 
used in this context in [Q. The initial and final states are defined by sets of complex 
variables a = {a\^}, b* = {b^}, respectively. A coherent state |ak) is an eigenstate of 
the annihilation operator ay^: ak|ak) = c^klo-k)- The relevance of coherent states for 
the scattering problem stems from the fact that upon differentiation they generate all 
momentum eigenstates (see eq.p!T| below). 

For simplicity, we illustrate the case of a single real massive scalar field with self- 
coupling g"^ in 3 + 1 dimensions. The results can be generalized to the case of gauge 
fields. First, we express the transition amplitude from an initial coherent state |a) 
at time Tj to a final coherent state |b*) at time Tf, in terms of the path integral by 
inserting a complete set of eigenstates of the field operator on the initial and final time 
slices (see figure 1) 



(b*|f/|a) = J d<j),d<f)f {h*\<j)f) {<f)f\U\<j),) {<j),\ai) 



(2) 
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Figure 1: Spacetime, with two slices. 



U is the evolution operator between time Tj and Tf. A "position" eigenstate of 
the field operator is denoted | 0) and 0ij = (j)(Tij). Then, from (Q), we obtain the 
S-matrix kernel in a compact form in terms of path integrals 



b*|S|a) = 5[b*,a] = lim 



d(j)f d(f)i e^i e^' 



D<p e 



iS[4>] 



(3) 



where S [0] is the action functional. The path integral appearing here is over fields 
obeying the boundary conditions 4>(Tij) = (pij. The functional Bf is 

Bf[h\cl^f] = -yd'kb*^b*_^e''-''''f -yd'kLO^<j>f0)M~k) (4) 
+ J d^k e"^'^f 6^ <Pf{-k) , 



tThe authors of § extremize momentum-space Greens functions. See Section 6. 



in terms of which the wave functional of the final coherent state is 

(b*|0;) = exp(i?^[b*,0;]) . (5) 
Similarly, the functional Bi is 

Bi [a, (Pi] = -i 1 d'k aka.ke-''"''^^ " ^ / "^'^ Wk</>i(fc)0i(-^) (6) 

is related to the wave functional of the initial coherent state 

(0i |a) = exp(5i [a, 0i] ) . (7) 
Here the 3-dimensional Fourier transform is defined 

Mk) = Jj^^''''^(T,j,x), (8) 

and is related to the residue of the 4-dimensional Fourier transform when the field 
reduces to a plane wave superposition as | Tj j | oo 

The kernel @ is a generating functional for S-matrix elements between any initial 
and final particle states, by functional differentiation with respect to arbitrary Ok 
and 6^. We now use this fact to construct a kernel for scattering from initial two particle 
states. We define an initial two particle (wave packet) state at t = Tj 

I p, -p) = J d^k anik) d^ j d'^k' aiik') aj^, | ) , (9) 

where d]^ is a creation operator, and q;_r,l(A;) are arbitrary smearing functions of fc, 
localized around some reference momenta p and —p respectively. The wave packets are 
normalized so that □ 

J d^k = 1 . (10) 

The regime relevent to the high energy multiparticle scattering problem is where, in 
the limit — > 0, the wave packet is localized around momenta \p\ ~ m/g"^ with a 
characteristic width, Ap ~ m. 

This state can be generated by functional differentiation of the coherent state |a) 
with respect to Ok 



\P, -P) = I d^k d^k' anik) aiik') 7^ | a) 

J OOk OClk' 



(11) 

a = 



So, differentiating under the functional integral, the S-matrix element between the 
two particle state ( |Tll) and any final state |{&k}) involves the following functional at 

t = T,, 



exp ( Bi [a, ( 



(12) 

a=0 



5ak 5ak 

2 .JUJ^, Uk) Hk') e-^K+-k')^» exp ( B^ [0, 0^] ) , 

after dropping a term which vanishes in the limit Tj —00. The last factor here is 
simply the normalization of the initial position eigenstate 

exp -i j d^k Uk Uk) U-k) ) . (13) 



■''We use the non-relativistically normalized commutator 0,^./ , 0^ — 5^ — k'^ 



We combine this with the smearing functions and finally obtain an S-matrix kernel 
for the scattering of two wave packets into arbitrary final states, 



S [b*, 2] = lim / d(j)f d(f)i an ■ (j), aL ■ 0i e^^ ['''^^] + f^' D<p e'^^'^'^ (14) 

where we have denoted the initial state ([TT|) by "2". Here we have used the following 
compact notation for the initial state factors 



a ■ 0i = I d'^k oi{k) (t)i{k) e-''^"^^ . (15) 



We now derive the boundary value problem obeyed by the classical field in the 
stationary phase approximation of (0). The stationary phase approximation of ordi- 
nary integrals suggests that the initial state factors be included in the extremization 
of the path integral i. These factors will certainly impact the stationary phase in 
the kinematic regime relevant to the high energy scattering problem, where the initial 
wave packets are peaked around \p\ ^ m/g'^, as expansions around the instanton indi- 
cate [^]. More generally though, we will show in Section 4 that this procedure provides 
a consistent weak coupling expansion at any fixed center of mass energy. 

We include initial state factors in a straightforward manner by first exponentiating 
the initial state factors into an "effective action", so that 

S[h\2] = lim f d(j)fd(f)iD(f) , (16) 

where the effective action F is 

r [0] = \naR■(|)iaL■(|)^ + 5, [0, 0,] + tS [(/>] + , (17) 

after dropping a term which vanishes as Tj — — oo. 

We can now derive the boundary value problem by varying the effective action. 
Varying the entire exponent F with respect to (f){x) for Ti < t < Tf gives the source-free 
equations of motion 

Varying the entire exponent with respect to (piik), gives 

^Mk) + ^kUk) = + ^1 e-^^' . (19) 

yaR- (pi OIL- (Pi j 

The first term on the left hand side comes from a surface term in the action S. The 
other terms come from variation of the wave functional at t = Tj. This boundary 
condition involves both the positive and negative frequency parts of the field, unlike 
the boundary condition which arises for an initial coherent state |a), which depends on 
the negative frequency component of only 0. 

The boundary condition (^) at the initial time slice is rather complicated. How- 
ever, it can be simplified since a real field may be written in the asymptotic region 
t = Tj — oo as a plane wave superposition 

Uk) = iu^.e-'-"^^ + wlkC^^'^^O . (20) 



^The Stirling approximation to the Gamma function is a classic example. 



Equation (|T9D then reduces to the requirement 



(21) 



(l + ^d^k anik) aiik) ) 



1/2 ' 



using the normahzation ([lOD . This solution is consistent with physical intuition, the 
classical field reducing to the initial particles at early times. The overlap of the left- 
and right- moving wave packets in the denominator is very small for narrow high energy 
wave packets. 

If the field is real, then the negative frequency part equals the complex conjugate 
of the positive frequency part, and the field (pi is determined as in ( pUD and (^TJ). The 
real initial condition can be integrated forward, and uniquely determines the final field 
0/. We say more about the case of complex stationary points of real fields in Section 5. 

It remains to be seen in what sense it is a good approximation to include the 
initial state prefactors in the stationary phase calculation. We will demonstrate in 
Section 4 that the stationary phase solution which results from this approach provides 
a controlled weak coupling expansion of the scattering amplitude. 

Finally, varying with respect to (f)f{k) gives 



We see that the precise scattering amplitude which is extremized is only determined 
at the end of the calculation, by the asymptotic form of the classical solution in the 
far future, Tf — ^ oo. For the initial conditions given by (0), the scattering amplitude 
corresponds to a transition from wave packet states to a final coherent state described 
by \b*) from (^2]). Note that neither a wave packet state i nor a coherent state are 
eigenstates of the Hamiltonian, so that neither have a definite energy. 

When the classical field satisfies the equation of motion (|18D and the boundary 
conditions (|T^ and (P^), there are no initial and final state corrections to leading order 
in the semiclassical expansion. This is not the case for the instanton, which satisfies 
the equation of motion but not the correct boundary conditions. Instead, there are 
linear terms in the fiuctuation about the instanton and these generate initial and final 
state corrections. 



The results of the previous section establish an important fact: 

The classical field obtained by starting with two incident wave packets and 
evolving them forward using (p!8| ) is the dominant contribution to some two 
particle scattering amplitude, in the semiclassical approximation. 

In the case of a purely Minkowskian classical solution, the result for the corresponding 
scattering amplitude will not contain an exponential suppression. Thus, if a real time 
classical trajectory can be found which connects wave packet initial conditions to an 
"interesting" final state, for example a multiparticle or fermion number violating state, 
our results imply that the corresponding scattering amplitude is unsuppressed. Classi- 
cal solutions which either exist in complex time or are themselves complex are relevant 
to classically forbidden transitions and will be discussed in Section 5. 

^Except in the limit of aji(k) ^ S^{k ~ p), where it reduces to a plane wave. 



-i^f{k) + uji,(f)f{k) = v/27j; 61k e'"*^^^ . (22) 
Again with a free field form as in (pO|), this boundary condition may be re-expressed as 




(23) 



SEARCH FOR SOLUTIONS 



The boundary conditions given by (|2^) and (^T]) are amenable to straightforward 
numerical integration. The boundary conditions specify 0(x,Tj) and (^(x, Tj), which 
are sufficient to construct 4>{x,t) for all subsequent t > Ti, given a discretization of 
the equations of motion ([18|) . In principle, it is possible to generate an infinite number 
of classical trajectories, each relevant to a particular initial wave packet scattering 
amplitude. 

Of course, it is not clear a priori that the final state which results from the classical 
evolution will be one which is interesting to the problem of multiparticle production 
or anomalous baryon number violation. However, the above observation summarizes 
the relevance for scattering amplitudes of computations by Rajagopal and Turok [^, 
and also Goldberg, Nash and Vaughn 0. Rajagopal and Turok studied the classical 
scattering of wave packets in the Abelian Higgs model, while Goldberg et al. studied 
the classical dynamics of 0^ theory. Neither group found that energy was readily 
transferred from high frequency to low frequency modes, a signal which would indicate 
the production of a final state with many (soft) particles. Therefore, the scattering 
amplitudes which are dominated by their classical trajectories are probably not relevant 
to multiparticle production or fermion number violation. 

On the other hand, the initial indications on the behavior of non- Abelian classical 
trajectories seem more promising. Some recent results on the classical behavior of pure 
SU{2) Yang-Mills theory by Gong, Matinyan, Miiller and Trayanov P| may provide 
insight into solutions of the boundary value problem presented above. These authors 
have considered the classical stability of a stationary mono-color wave in Yang-Mills 
theory (in Aq = gauge): 

A^{x,t) = 6c3 A cos/cq^; cosciJot, (24) 

where c is a color and i is a spatial index 0. Small amplitude variations of the field in 
directions of different color are found to lead to an instablility with long wavelength. 

The Hamiltonian evolution of a perturbed standing wave is equivalent to the evo- 
lution of the initial conditions relevant to the scattering of two perturbed plane waves 
in SU{2) gauge theory. The discovery of an instability of the initial configuration to 
decay into long wavelength modes implies the existence of a classical trajectory con- 
necting initial high energy plane waves to long wavelength modes in the final state. It 
therefore suggests that the corresponding 2 many gluon scattering amphtude may 
be unsuppressed ! The relevance of weak coupling calculations to pure gauge theory 
is unclear, due to the asymptotic freedom of the theory. However, it is possible that 
arguments similar to those from deep inelastic scattering may be used to justify the 
semiclassical result. In the limit that all of the energy scales (including the energies of 
the many "soft" outgoing gluons) are large compared to the intrinsic mass scale of the 
theory h.su{2)i it seems plausible that the relevant running coupling constant is small. 

It is also important to determine whether the instability persists for wave packet 
initial conditions, as the plane wave limit is never achieved at an actual accelerator. 
Since wave packets are localized in space, two packets have only a finite amount of time 
to interact before passing completely through each other. This is in contrast to plane 
waves, which exist everywhere in spacetime. Since the plane waves of Miiller et al. |Q 
have finite, nonzero amplitude, they also contain the equivalent of an infinite number 
of particles. It is quite possible that the behavior of wave packets representing a finite 
number of particles is very different. 

Another class of trajectories in pure SU{2) Yang-Mills theory are those found 
by Farhi et al. 0. These solutions in the spherical ansatz correspond to spherical 
shells of glue that originate at infinity, collapse inward, and bounce back to infinity, 
leaving behind some fractional topological charge. There are some open questions 
concerning such solutions - in particular, their stability to perturbations which are 

II NB - A travelling wave is stable in Yang-Mills theory 



outside the spherical ansatz, and their relevance to fermion number violation. However, 
in the absence of fermions, they can be interpreted within our formalism as extremizing 
scattering amplitudes between spherically symmetric initial and final quantum states. 

Of course, our eventual goal is to address the electroweak theory, where the Higgs 
mechanism cuts off the infrared growth of the coupling constant. It would be ex- 
tremely interesting if simulations of wave packet scattering could be performed in a 
spontaneously broken SU (2) gauge theory. 

Clearly, there is much numerical work to be done to address these issues. 



THE EXPANSION 



Let us now construct the perturbative expansion around the classical solution of 
the EVP. We will demonstrate that corrections to the leading contribution are system- 
atically suppressed in powers of the small coupling constant for any fixed energy. The 
resulting small coupling expansion differs from the usual one (e.g.- instanton) and the 
subtleties involved will be discussed at the end of this section. 

As usual, we proceed by expanding the quantum field around a classical back- 
ground field (j)'^, which is a solution to the classical boundary value problem derived in 
Section 2. 

0(x) = (j)''{x;n) + u{x) . (25) 

The classical solution depends in general on a collection of collective coordinates de- 
noted here by /x, corresponding to the invariances of the scattering amplitude 0, 0]. We 
expect invariances of translation and scale size to be broken by the choice of (yR^iik)', 
however, /i may include gauge coordinates if we are working in a gauge theory. Equa- 
tion (p5D involves an expansion on the initial and final time slices, which we can express 
in momentum space as 

<PiAk) = <Plf{k;f^) + ^^J{k) for t = Tij, (26) 



where (1)1 j{k) is known explicitly for any given initial and final wave packets, from ( poD 
and (|2|)' 

Note that the classical solution is of order 0{g^), since it is non-trivial in the limit 
as (7 — 0. In fact, when g vanishes, the background field is simply the initial condition 
propagated forward in time with the free Hamiltonian, 

<P\t,k) = («ke-'-^* + <ke'"^*) , (27) 

where is as in (PT]). In this limit, the two wave packets simply pass through each 
other. For non- vanishing g, the classical field (p'^ will have some complicated dependence 
on g, reflecting the non-linearity of the theory and the external boundary conditions. 

The kernel of the S-matrix (0) is composed of three functional integrals which we 
write as 

S [b*, 2] = / d^f d<P, an ■<P,aL- (t>i e""^'"*'*^' + ""^'"'^^^ ^ (28) 
J j<j>i 

where the limit |Tjj| — oo is understood. Now we formally expand each of the terms 
in (p8[) using (|^). The action becomes 

S[cP^ + u] = S[<P'] + (d^xu{x)^\x) + S2[u] + S,nt[i^] , (29) 



after integrating by parts and retaining the surface terms in the time direction, and 
using the equations of motion (|18D . All of the non-quadratic dependence is contained 



in the interaction terms Smt — 0{v^,v^) which are suppressed by powers of g^. Here 
5*2 is the part of the action quadratic in the fluctuation field, which can be expressed 

S2[y] = -I f d^x u{x) z>(x) * - ]- I d^x v{x) A [x] u{x) , (30) 
2 J Ti 2 J 

in terms of the operator A of quadratic fluctuations in the 0^ background, 

6S 



5 (x — y) A\x] = - , , , ^ , , , 
Expanding the boundary functionals gives 



5\x-y){dl + V"i<P^)} . (31) 



-i / d-'k Uf{k) <p%-k) + Bf 



after using the boundary condition (^2]) , and 

I [ d'k Uiik) ^'r{-k) + Bi[0,(j)^] + Si[0,Z/i] - 



(32) 



(33) 



after using the boundary condition ([T9|) . The 0jj terms in (^) and (|33D cancel with 
the boundary terms in the action (pQl). Note the appearance of additional terms linear 
in i^i in (|33D , which arise because we are expanding around a stationary point of the 
effective action (0), instead of the action. Note also that the terms independent of Vij 
are the overlaps of the classical field with the initial and final states, respectively. 
Now, substituting (|2|), (|32|) and (|3]) in (|28|), we obtain 



^[b*,2] 



(34) 



where the quadratic integral over fluctuations is contained in 

Zm = (35) 



after using ■ 0^ = 1 in the initial state factors and in the exponent. An integration 
over collective coordinates has been factored out of (0) and the remaining functional 
integral ( P^D involves only fluctuations orthogonal to any zero modes of the classical 
background. 

The 0{g^) corrections result from the interactions in Sm- Expanding exp{iSint[i^]) 
inside ( pSD yields the leading term plus corrections which have the form of vacuum to 
vacuum loops. The vertices used to construct these loops each carry a suppression of 
g^. The propagators inside these loops are the usual Feynman propagators evaluated 
in the classical background (j)'^ ta. In the instanton case 0, |Tl| , dangerous initial state 
corrections arose from the residues of propagators in the instanton background, which 
displayed an energy dependence ~ g'^s, where s is the center of mass energy of the 
collision and could be parametrically of order 1/g^. Note that it is only the residue 
of the instanton propagator that is ill-behaved and not the propagator itself. Since 
the corrections in our case are proportional to loop integrals rather than the residues 
of individual propagators 0, |TI|, we expect them to remain subleading even at high 
energies. 



**It is worth remembering that since we are working directly in Minkowski space (instead of be- 
ginning in EucUdean space and analytically continuing back), the Feynman boundary conditions on 
Greens functions are only obtained through some form of regularization, such as an ie prescription. 



The Gaussian path integral appearing in ( ^51) can now be done exactly. First factor 
out the quadratic prefactor by defining a more general functional integral, 

Z [j] = I dvf du, Du e^/P''^/] + ^'["''^^l + t'^] - I'''' , (36) 

in terms of which (^) may be written 



IJ=J* 



where we will set the arbitrary current j{k) equal to 



= iMk) + aL{k)) e-^-k^» , (38) 

after the functional differentiation. The stationary point of ( |5BD is a "classical" fiuctu- 
ation field z/^ satisfying 

A [x] z/^ (x) = , (39) 

with boundary conditions 

iu^{k) + uji^u-ik) = j{k) and -iu'f{k) + uJi,u%k) = 0. (40) 

The solution of this linear homogeneous boundary value problem can be expressed in 
terms of a Greens function 



iy{x) = J d'qGix,q;t)j{q) , (41) 

where G obeys 

A[x] G{x,x';t) = , (42) 
with boundary conditions in terms of its Fourier transform 

tG(k,q;TA + uJi,G(k,q;Ti) = 6' (k - q) (43) 



and 

~tG{k,q;Tf) + U},G(k,q-Tf) = 0. (44) 

Using the above expressions, we obtain a compact expression for ( PBD entirely in 
terms of this Greens function, evaluated at t = Tj, and the known initial state source. 

Z [j] = det-i/2 [A] exp Q | d'k d'q j{k) G {k, q; T,) j{q) ) . (45) 

The exponent appearing here is due entirely to the boundary functional Bi in (|36|) . All 
other terms vanish on the stationary phase solution (|^). With this result, we obtain 
the Gaussian integral (^) by functional differentiation and setting j = 

Z[(j)'] = det-^/2[A] . (46) 

Several terms in (^) have vanished in the limit of Tj — »• oo, most notably initial 
state factors. This is because both j*(/i;) and our inner product (|T3|) contain rapidly 
oscillating factors e"*"^*^' . This rapid oscillation allows the application of the Riemann- 
Lebesgue lemma in the limit Tj — * — oo, which guarantees that for any function f{k) 
whose Fourier transform exists, 

lim / dk f{k) e-''^"^^ = 0. (47) 



Here the functional determinant det^^^^ [A] is actually the product of three deter- 
minants: the standard 4-dimensional determinant of the operator A satisfying Feynman 
boundary conditions, and two 3-dimensional determinants representing the edge fluctu- 
ations v{x,Tij). The latter determinants can be expressed in terms of a homogeneous 
Greens function similar to G described above. 



We have demonstrated that the kernel of the S-matrix p8| ) is obtained to leading 
order in from the action of the classical solution and the determinant of quadratic 
fluctuations around it. 

S[h\2] = /rf/xe^^'[^''^/l + ^^[°'<^»^l + '^l^^Uet-i/2[A] \l + 0{g^)] . (48) 



The classical factors here are known functions, depending only on the initial and final 
boundary values of the classical field, ( pOD and (pSf). Thus, we may further reduce ( ^S]) 
to 

S[h\2\ = e-' Jdfie^/' + '^^'^'Uet-'/'[A] [l + 0{g^)] . (49) 

where N = J d^k 6^6^ is the average number of particles in the final coherent state, 
which depends implicitly on the classical field. The origin of this factor is the use of 
coherent states which are not normalized. The standard normalization for coherent 
states which we use here is 

(b*|b*)=e^. (50) 

The correctly normalized S-matrix element would not include this factor, leaving the 
classical action and the determinant. 

There are no additional tree- level radiative corrections to (^), as there are in 
an expansion around a stationary point which does not obey the correct boundary 
conditions (c.f.- the instanton 0, ^). This is a direct result of the fact that our 
classical solution obeys both the correct equation of motion and the correct boundary 
conditions. 

It is worth noting that our semiclassical expansion differs from the usual one in a 
crucial way. Our classical background depends on the coupling constant itself. In the 
usual instanton calculation, the field equations and boundary conditions can be made 
independent of g by rescaling — >• = gcp. In our case, since the boundary conditions 
are fixed by the wave packet shapes aR^iik), the equations of motion for contain a 
dependence on g which cannot be scaled away. It is clear that a change in g will lead 
to a different nonlinear evolution of and hence a different classical trajectory. In 
the limit of strictly zero coupling constant, the evolution of the initial wave packet is 
clearly trivial, and we get only a contribution to the "diagonal" part of the S-matrix 
(i.e.- no scattering). Therefore, in order to obtain a nontrivial result, we have to 
take g large enough to find an interesting classical trajectory, and yet small enough 
to justify the semiclassical g"^ expansion derived above. We have not demonstrated 
that the two requirements can be satisfied simultaneously in any theory. This can only 
be checked a posteriori by computing the size of the subleading 0{g'^) corrections in 
the background of a proposed trajectory. Of course, our experience with perturbation 
theory leads us to believe that for g"^ <^ 1 the corrections will be small. Thus, a 
classical trajectory found using weakly coupled equations of motion is likely to provide 
a controlled approximation. 



COMPLEX TRAJECTORIES 



In the preceding discussion, we have focused on classical trajectories dependent 
on a real time variable. These are relevant to classically allowed amplitudes, such as 
those for multiparticle scattering or transitions at energies well above an energy barrier 
in configuration space. As we noted previously, the existence of a real time trajectory 



implies an unsuppressed scattering amplitude for the corresponding transition. How- 
ever, there are many instances in which classically disallowed transitions are of interest. 
In particular, baryon number violation in the electroweak theory at energies less than 
the sphaleron energy Eg ~ M^/aw ~ 10 TeV is classically forbidden and therefore 
cannot be described by the same type of real time trajectories. The application of our 
formalism to such cases therefore requires some generalization, which we will describe 
in this section. 

Consider the initial value problem given by ([TSD and ([T9|). Since both the initial 
field values 4>i{x) and their time derivatives (f)i{x) are given, the classical field has a well 
defined energy. Since real time trajectories conserve energy, it is clear that a suitable 
trajectory will not exist for some choices of initial and final conditions. In particular, in 
the electroweak theory for < Eg, there are no real time trajectories which connect 
initial and final conditions with different topological charge. 

In order to find such trajectories, previous authors [|T^, |T3|, [T^ have proposed com- 



plex time paths, which involve alternately Minkowskian, Euclidean and Minkowskian 
paths joined together at two times, say and iT as in figure 2. In this approach, t = 
and t = iT must be turning points %) = for all x) of the field equations with 



given initial conditions. As is familiar from one dimensional quantum mechanics [ITS 
complex time contours describe an incident configuration which reaches a turning point, 
and then tunnels under the barrier in Euclidean time, returning to real time at the sec- 
ond turning point upon reaching the other side of the barrier. The field (f){x) can 
be taken to be real on the entire contour because of the existence of the two turning 
points tl. The action is then purely imaginary on the Euclidean part of the contour 
and yields the familiar exponential tunneling suppression. 

\{f\^)\' - e~''-' . (51) 

As a matter of practice, one might try to find a suitable complex time trajectory 
by integrating our initial conditions (^) forward in real time in hopes that a turning 
point is encountered. Then, the time variable should be Wick rotated t 1x4 and the 
turning point configuration integrated forward in Euclidean time until another turning 
point is encountered. Finally, the contour should be rotated back to real time and 
the second turning point integrated forward until it achieves its asymptotic plane wave 
state. The final state \b*) can be read off from (p3D. 

Im t 



iT 



Re t 



Figure 2: Complex time contour 

However, in the most general case, there is no guarantee of encountering turning 
points. Here we propose a more general method for extremizing S-matrix elements 
for which there is not a corresponding real, Minkowskian trajectory. Rather than 
complexifying time, we instead search for a general complex saddlepoint configuration 



satisfying ([T8|) and (0). Our motivation stems from the usual method of steepest 



t^A fact guaranteed in one dimensional QM, but not in higher dimensions or in field theory IT 



descents applied to a real integral of the form 



I{h) 



(52) 



In general, the asymptotic series in h for this integral can be obtained by expanding 
about a saddlepoint of the function /(x) : f'{x)\x* = 0. However, the saddlepoint x* is 
often complex, and to apply the method, one must first deform the original real integral 
into the complex plane (see figure 3). 

In our case, we first think of the Feynman Path Integral as the product of a large 
number of regular integrals by discretizing spacetime: 



n d'^x 



(53) 



For a real field 0(x) each integral in the product is a real integral like (0). We generalize 
the method of steepest descents by allowing the variable of each integration, to 
become complex. There should be no obstruction to this generalization as the function 
S{^x,^x±i) is analytic in the variables This procedure leads us to search for a 
complex trajectory {^x} (or in the continuum, (f>'^{x)) satisfying (p]8|), ( |T9| ) and (|2^). 
The integral over small fiuctuations about (j)'^{x) must be performed along the path of 
steepest descent in configuration space. In other words, along small z/(x) fiuctuations 
which keep the imaginary part of iS [(f)'^{x) + z^(a;)] constant. 



Im X 




Re X 



Figure 3: Contour in complex plane, 

deformed through saddlepoint x* 



For a complex field, there is no longer a relation between the positive and negative 
energy Fourier components. Without loss of generality, we can write 



2Ur 



Up e'''^'^^' + Vp e'^"'^' ) 



(54) 



where Up , Vp are independent complex functions. The initial boundary condition (|T9|) 



is satisfied by 



ttij(p) 



Up = '^^'''L + "^^^l , (55) 

/ d^k aR{k) t>k / d^k aiik) fk 

where is arbitrary. The initial configurations which satisfy ([5^), ([55| ) have negative 
frequency modes which "resemble" the wave packet ai + an (up to complex multiplica- 
tive factors) and positive frequency modes which are arbitrary. 

It is then clear that there are an infinite number of complex trajectories which can 
result from evolving ( 15^) using ([T8|). The final quantum state |b*) which results from 



a given trajectory depends only on the positive frequency component of the trajectory 
in the far future (see Eq. (22)). The multiphcity of trajectories allows a very large 
set of scattering amplitudes to be addressed within our formalism. In fact, it seems 
in general there always exists a complex trajectory which satisfies ([TSD , (p!9|) and ( p2D 
simultaneously for any a{k) and b*. This is because the initial and final conditions 
(0) and (|22| ) each correspond to conditions on a combination of (j){x)\T, and 0(x)|Ti, 
and 0(x)|t^ and <P{x)\tj:, respectively (mixed Dirichlet-Neumann boundary conditions). 
Since the equations of motion (^) are second order in time, this mixed set of boundary 
conditions should be enough to specify a unique solution. 

Since the trajectories found here are complex, the action «S'[0'^] may contain a real 
suppression factor. This is to be expected, as some of the trajectories will correspond to 
classically disallowed tunnelling transitions. Thus, although potentially any S-matrix 
element can be approximated in this way, the result for particular trajectories may be 
exponentially small. 



RELATION TO PREVIOUS WORK 



Here we comment on the relationship between our approach and some previous 
work by other authors on high energy multi-particle and anomalous processes in elec- 
troweak theory. First, we will recall the expansion around the constrained instanton in 
electroweak theory This approach neglects effects of initial and final states on 

the saddlepoint and results in a low energy {E <^ E^, c:=l M^/aw) approximation to the 
total two-particle cross-section. Then, we remark on a previous approach to account 
for the impact of initial and final states on the saddlepoint and comment on the 
formal resemblance which it bears to our approach. 

A few of the properties of the electroweak instanton will make clear the relation- 
ship to our work: {i) The electroweak instanton is a solution of the Euclidean equations 
of motion. So, it will be necessary to rotate our real time formalism to Euclidean space 
in the analysis above, t — > ix^. We would then be considering a saddlepoint approxima- 
tion to an on-shell truncated Euclidean Greens function, (m) The electroweak instanton 
has finite Euclidean action, S = Sir'^/g'^, and therefore satisfies vacuum boundary con- 
ditions. So, it is a suitable saddlepoint point for a Greens function with external fields 
only insofar as the effect of initial and final states on the saddlepoint are neglected 
(dropping the right hand side of (|19D and (^)). 

Since the instanton does not obey the correct boundary conditions to be a sad- 
dlepoint of a scattering amplitude @, linear terms in the fluctuation expansion do not 
cancel, and the expansion of the scattering amplitude entails corrections which are for- 
mally large 0{1/ g"^) in the exponent. Fortunately, these corrections are under control 
and calculable for sufficiently small energies, E <^ E^, The effect of the final states 
can be taken into account, in a perturbative expansion in powers oi x = E / E^ <^ 1 [Q. 
For instance, the total inclusive cross section in the one-instanton sector is given by the 
well-known result 

<yUE) ~ (-1 + <^i^"" + C2X^ + O (x'/'))^ , (56) 

in the limit E/E^, = fixed and g 0. Here ~ implies that only the exponential 
behavior of the cross-section is shown. The first term here is just twice the instanton 
action, the 'tHooft suppression for vacuum tunnelling, while the next term indicates the 
exponential growth of the cross-section at low energies The higher order terms are 
determined from all tree-level corrections to the many soft {E ~ My^) final state parti- 
cles, and the first few coefficients q are known. Thus, provided the energy is sufficiently 
low E E^, tree-level corrections to the final state can be described semidassically. 

It has not been similarly demonstrated that the effect of initial state corrections, 
involving many loops, can be described semidassically. However, some results indicate 



that these corrections may run counter to naive intuition and the initial two particle 



state may be described semiclassically as well [|TT], 0, 3. As we noted in Section 4, 
initial state corrections in the instanton expansion appear as factors of the residue 
of the propagator in the instanton background. The residue of the propagator has 
hard high energy behavior ( ~ g'^s ) where s is the center of mass energy |]1T|, so that 
initial state corrections can contribute at the semidassical level when s is of order 
Mueller [|ll| has shown that they contribute first at order (E/E^)^'^^^ to (|56D 



By contrast, in our formalism the classical background field is constructed to solve 
the correct multi-particle boundary value problem derived in Section 2. So, there are 
no additional tree-level corrections in the limit E = fixed while g ^ 0. All corrections 
are formally suppressed by powers of g, as shown in Section 4. The size of subleading 
corrections to our result depend on the properties of the Feynman propagator in the 
background of our proposed classical trajectory, (p'^. It would therefore be very interest- 
ing to understand the high energy behavior of the propagator in our background field. 
Also, for the purposes of direct comparison to (pGf), it would be interesting to obtain 
its behavior in the same limit, where E/E^ = fixed while g 0. We have reason to 
expect that the propagator behaves less severely than the instanton propagator in this 
limit, due to the lack of translational zero modes of the background field 

A previous approach to account for the impact of initial and final states on the 
saddlepoint P] bears a resemblance to ours but results in a different classical equa- 
tion. We discuss it here for the purposes of comparison. These authors extremized 
Minkowskian n-point Greens functions of electroweak gauge fields 

G^") (pi,p2, . . .) = f'^^ ^^'"Hpi) ■ ■ ■ e'^t^l/^' . (57) 



It is not clear to us whether the Greens function is the proper quantity to extremize 
to yield information about the corresponding scattering amplitude (and ultimately the 
cross-section). It may be that extremizing the entire Greens function, rather than its 
LSZ projection, is too strict a requirement. 

An equation of motion for the saddlepoint of (pTf ) can be derived in a manner with 
formal similarity to our own. To do so, first express the n fields as an exponential, 
using A = explnA, and then extremize the exponent. This yields 

■ y 5^^" 5'''' -A—- . (58) 

' " /I Ll'i Qi'i ( fy-t . 1 



This is a non-linear integro-differential field equation, depending on both the field A{x) 
and its Fourier transform A{p). Not much is known about equations of this type, and 
despite much effort (|58| ) defies solution [jl5| . It appears more difficult to solve in practice 
than the boundary value problem we have derived (0), (^) and (|22|) . 

Before concluding I would like to point out a difference between studying scattering 
via the Minkowskian and Euclidean functional integral. The Euclidean counterparts 
to our solutions (i.e. their analytic continuations to imaginary time t ix^) are very 
badly behaved asymptotically. That is, the 'wrong frequency' components required by 
the prescence of particles in the initial and final states (see equations (|1^) and ( p2D ) 
imply an exponential blow-up of the solution at large imaginary times. In other words, 
our BC's are incompatible with the usual Feynman BC's (only positive frequencies in 
the future, and negative in the past) which arise automatically from requiring finite 
action in Euclidean space. 

SUMMARY 



In this lecture I have outlined a connection between certain classical solutions 
and corresponding quantum S-matrix elements in a weakly coupled field theory. The 



connection is made via a semiclassical approximation which appears to be controlled at 
small coupling, regardless of the scattering energy. This is in contrast to methods which 
involve expansions about instanton-like trajectories (i.e.- satisfying vacuum boundary 
conditions), which have been shown to break down at nonperturbative energy scales 
E ~ m/g"^. Real Minkowskian trajectories can be found by the straightforward time 
integration of a well-defined initial value problem which is determined by the inital 
quantum state. Complex Minkowskian trajectories require the implementation of mixed 
boundary conditions at both initial and final times. This procedure, while more difficult 
in practice, can be guaranteed in principle to describe any wave packet to coherent state 
transition specified by the complex functions a{k),bl. The formalism described here 
can thus be used to compute intrinsically nonperturbative scattering amplitudes in a 
variety of field theories. 

Some examples discussed here are the problem of baryon number violation in 
the electroweak theory and multiparticle scattering in gauge and scalar theories. The 
discovery of nontrivial or "interesting" real classical trajectories in these theories can 
now be related to unsuppressed scattering amplitues which are potentially observable 
experimentally. Alternatively, nontrivial complex trajectories, of which there are an 
infinite number, allow the calculation of a much larger class of scattering amplitudes, 
including some that are classically forbidden and therefore involve tunneling. We hope 
that our results will stimulate future numerical work in this area, particularly the search 
for nontrivial trajectories. 
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